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Abstract 

A product difference equation is proved and used for derivation by elementary 
methods of four combinatorial identities, eight combinatorial identities involv- 
ing generalized harmonic numbers and three combinatorial identities involving 
traditional harmonic numbers. For the binomial coefficients the definition with 
gamma functions is used, thus also allowing non-integer arguments in the iden- 
tities. The generalized harmonic numbers in this case are harmonic numbers 
with a complex offset, where the traditional harmonic numbers are a special 
case with offset zero. 
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1 Introduction 

Binomial coefficients and their combination with harmonic numbers occur frequently 
in applied mathematics [TH [T31 [03 [35] . Many ways to find and prove these identities 
can be found in literature [5[ [1(3 H3 [23 [HJ [231 113 • In this paper a product difference 
equation is used for derivation of a number of such, often new, identities. 

2 A Product Difference Equation 

Theorem 2.1. Let n be a nonnegative integer, and let x and y be complex variables 
and let the set {zk} be n complex variables and the set {uik} be n complex variables. 
Then there is the following product difference equation: 



k=l fc=l 
n k—1 n \ ' / 



ll(x-z k )^-l[(y-z k ) 

k=l k=l 
n 

= ]T[(.t - Zh r* -( y - z k ) wk ] i[(x Zl r U(y- zi) 

k=l 1=1 l=k+l 

where a sum of zero terms is zero and a product of zero terms is one. 
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Proof. Let a k = (x - z k ) Wk and /3 k = (y - z k ) Wk . Then ^27TJ> becomes: 

n n n k—1 n 

flak- JJ^ fc = 53(^-/3*) JJ««- II A 

fe=l fc=i fe=i z=i Z=fc+1 ^ 

n fe n k—1 n 

=Edi««- n A-n^-iiA) 

fe=l ;=i i=fe+i z=i z=fc 

On the right side the first product for each k cancels the second product for each 
k + 1, only leaving the first product for k — n and the second product for k = 1. 
These two remaining products are the two products on the left side. □ 

In the case that all z k = and all w k = 1, (|2.ip reduces to the following power 
difference equation [19]: 

n 

x n -y n = (x-y)Y / x k - 1 y n ~ k (2-3) 
fe=i 

3 Combinatorial Identities 

The definition of the binomial coefficient in terms of gamma functions for complex x, 

y is PEJEI1EE2]: 

'*\ = r ^ + 1 ) (3 1) 

K yJ T(y + l)T(x -y + 1) V ' ; 

For nonnegative integer n and integer k this reduces to [U El [T2l [13] : 



if < k < n 



= {k\{n-k)\ ~ ~ (3.2) 

otherwise 



From the recurrence formula for the gamma function for complex s [T] [5] [TTJ [TH1 [M] : 

r(s + i) = s r(s) (3.3) 

follows for integer a < b and complex s: 

Combination of (12 . 1|) with the {z^} consecutive integers and all w k — w and (|3.1[) and 
one of the right side expressions of (13.4[) for a; and for ?/ in (|2.ip yields the following 
four combinatorial identities. 
For integer a < b and complex w, x, y: 



f x + k 



J2[(x + k + l) w -(x-y + kr] 



fc= 



V 



(y + i) 1 



x + b\ w fx + a 
y + l) _ U + 1 



(3.5) 
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6-1 



J2[(y-x + k) w -(y + k) w ](-i) wk 



k—a 



u> b 



x - 1 
y + 6-1 



(-1)« 



a; 

2/ + fe. 

x - 1 
y + a- 1 



6-1 



w k 



j2K-y-k-ir-(x-y-k + ir}(-i) 

x + l 



k—a 



X 

y + k 



(x + iy 



(-i) 



wb 



x+l 
y + b 



y + a 



6-1 



J2[(x-kr-(y-k + iy 



x\ /y> 



(y + i) 1 



x\ (y + 1 



x \ /y + 1 



For to = 1 these identities reduce to 

6-1 



6-1 



6-1 



X 

y + fc 



£\ / y N 



E 

k—a 



k \k 



E 



X 

y + k 



_ x + 1 
~ i + 2 

y+1 
x-y- 1 



x + k\ fx + b 



x + a 



\y + y \y + 1 



-i) a 



x-l 
y + a-1 



(-1)« 



X+l 

y + a 



(-1)' 



(-l) f 



x-l 
y + 6-1 

x + l 
y + 6 



x\ /y + 1 



A /y+l 
a J V a 



which are equivalent to known identities [?]■ 
For an example of iy = 2, let be given [7j IT3"]: 



E 

fe=0 



/2n 
\ n 



Then ([3"1>]I with «; = 2, a = 0, b = n+l, x = n, y = and (f3"T5")) yields: 



E* 

fe=0 



n \ n 



2n 



2 V n 



For an example of u> = 3, let a formula from A.C. Dixon be given 



In , s 3 



E(-d a 



= (-iy 



2ra\ /3n 
n \ n 



Then with w = 3, a = 0, b = 2n + 1, x = 2n, y = and (j3~T5|) yields: 

B-l)^(2~- t )f 2 :) a = (-ir^f 2 „")( 3 ;) (3.16) 



4 Combinatorial Identities with Harmonic Numbers 

The definition of the harmonic numbers for nonnegative integer n is |13[ 115] : 

- 1 

H n =J2z (4-1) 



k 
fc=i 



from which follows that Hq = 0. The definition of the generalized harmonic numbers 
for nonnegative integer n, complex order m and complex offset c is [141 115] : 



1 

= Y 7 7T- (4.2) 

^(c + k) m v ' 

k—l 

from which follows that H^q = 0. There is the following symmetry formula: 

i?M = (-l) m H (n ! ] ^ , „ (4.3) 

For order m = 1 and nonnegative integer offset c — k these numbers can be expressed 
in the traditional harmonic numbers of (|4.1[) : 

= H k+n - H k (4.4) 

from which follows: 

H$l = H n (4.5) 

Identities involving these numbers can be obtained by differentiation of the product 
(|3~4"|) and using 

± f[ (x kr = V« n(x - fc)» = - k) w (4.6) 

and similarly for y. These identities therefore are in pairs: one from differentiation to x 
and one from differentiation to y, in each case using one of the right side expressions in 
(|4.6p . This results in the following eight combinatorial identities involving generalized 
harmonic numbers. 

For integer a < b and complex w, x, y: 



X> + k + l) w H^ k _ a+1 -(x-y + kYH^ k _ a ] f:r + ' 

k—a 



(4.7) 
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/ x + k 



k—a 



y 



y _|_ I J y-x-b,b-a 

6-1 



k—a 



•(-i) 



6-1 



wfe /" ^ A _ rr wf_i \wb ( X 1 \ (i) 



Et(l/ - ^ + fc)"'^ ) -6 > 6- fe -l - + fc) W ^-6,6- fe ] 
fe= a 



•(-lW M =-a;" , (-ir a f * 1 ) 
+ \y + a-l> 

6-1 



y — 6,6— a 

]T[(-y - k l) w H^ k _ a+1 - {x -y-k + l) w H^ k _ a ] 



k—a 

■ (-iy 

6-1 



w / _j_ -1 \ — «> 



A;— a 

■ <- i »" c + *r ■ _(i + ir( - ir g + o 

6—1 / \ w / \ ~~ 10 

5> - fcr^-l-M-a+i -(»-*+ ir^-U,*-.] (I) f I) 

k—a ^ ' ^ ' 

6—1 / \w / 

Y> fc ) tu <Vi,6- fe -i -(»-*+ ii^S+u-j Q ( 



\ w / 1 1 

X \ / 7/ + 1 



For io = 1 these identities reduce to: 

6-1 



V ( X + k \H {1) 

M y / - 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



6-1 

E 

k—a 



x + k 

y 



H. 



(i) 

y — x — b,b— k — 1 



X + ' (H {1) 



-) 



1 /x + b 



6-1 



k—a 



(-i) 



y + k 

x - 1 
y + 6 



(i) 



y— x+a— l,fe— a 



1 W (1) + I )+ (zDV 

-l) K !/-x+a-l,6-a ^ ^ ^ + 



-6,6-fc-l 



r(i) 



(-i)V 



x - 1 



x \y + 6 — 1 



6-1 

/c— a 
X + 1 



(1) 

y+a.A; — a 



x + 2 

6-1 

£(-!)■ 

/c— a 
X + 1 



(-1) 



6+1 



1 _ (-l) a ^x + 1\ 1 



( X 

\y + k 



( x + 1 Y\ H w 

\y + bj y v+a,b-a X + 2 J x + 2 \y + a 
1 



H 



(i) 

x — y — b-\-l,b—k—l 



x + 2 



(-1) 



6-1 

E 

fc— a 



(i) 



l x-y-b+lM-a x + 2 > x + 2\y + b) 



y + 1 



k)\k) a-x-l.k-a x -y-l 



x\ (y + 1 



K^a-x-lfi-a ^ 



1 



-) 



1 



6-1 

E 

k—a 



x\ y 
k \k 



x\/y + l 
a) \ a 



H, 



(i) 



x — y — 1 x — y — 1 \a 
y + 1 



x\ ( y + 1 

a 



y-b+hb-k-i y - x + i 
-l 



(1) 

y— 6+1,6— a 



i 



1 /x\ (y + 1 



y — x + 1 y — x + l\b 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



(4.22) 



Changing all k into a+6 — — lin the summation term reverses the order of summation 
and yields valid identities. When the offset of the generalized harmonic numbers in 
these identities is zero, these identities yield one traditional harmonic number identity 
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and three combinatorial identities involving traditional harmonic numbers. 
For nonnegative integer n and complex m, w: 



[(* + l)™ - k w ] H k = (n + l) w H n+1 - Hg£$ (4.23) 

k=0 



Y,( k )H k =( n + 1 \ Hn+1 -^) + ^-( ° ) (4.24) 
^\m) \m+ir + m+l' m + 1 V m + 1 ' 



fe=0 



k=0 



-1 



n — k J \k J n — m+l 



k=0 



to \ m+l 



D-d*" a- 



to + 2 



1 / / TO , 

ffn+l + — T "I 

n — to + 1 V n + 1 



to + 1\ X ,„ 1 



(4.25) 



(-1)" ' (ffn+l- 



n + 1 / to + 2 to + 2 



(4.26) 

The identity (|4.23l) is an example of harmonic number identities with polynomials or 
rational functions in k [20l [22] . The identity (|4.24l) for nonnegative integer to was 
already listed and proved in literature [TU [13] . More complicated identities involving 
harmonic numbers with higher orders may be derived from the given identities by 
using (|4.3[) and the following relations [31 U] [5] [7] . 
For nonnegative integer n and complex to, w, x, y: 



dx 



H^ n = -mHg+$ (4.27) 



d(x + y\ w _ {1) fx + y 



dx\ n x+y-n^n 



II 



= wH^'J ^ a (4.28) 
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